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I. INTRODUCTION 



The lossy Bosonic channel provides a quantum model 
for optical communication systems that rely on fiber 
or free-space propagation. For the pure-loss case, in 
which the quantum noise accompanying the loss is the 
minimum permitted by quantum mechanics, the clas- 
sical information-carrying capacity of this channel has 
been derived, and shown to be achievable with single- 
use coherent-state encoding 0. For the more general 
thermal-noise channel, in which the environment injects 
an isotropic Gaussian noise, the Holevo information of 
single-use coherent-state encoding is a lower bound on 
the channel capacity that is tight in the limits of low and 
high noise levels 13 Moreover, if a recent conjecture 
concerning the thermal-noise channel's minimum output 
entropy is correct, then single- use coherent-state encod- 
ing is capacity achieving £|. 




c = jrj a + sjl-rjb 



FIG. 1: Two-user, single-mode, optical multiple-access chan- 
nel. Transmitters Alice and Bob have access to input modes 
a and b, respectively. Charlie receives the output mode 
c = y/rja + y/l — nb. 



II. COHERENT-STATE MAC 



To date there has been almost nothing reported 
about the classical information-carrying capacity region 
of multiple-access Bosonic channels, i.e., Bosonic chan- 
nels in which two or more senders communicate to a 
common receiver over a shared propagation medium. In 
this paper we derive single-mode and wideband capac- 
ity results for such channels |]|. First, we show that 
single-use coherent-state encoding with joint measure- 
ments over entire codewords achieves the sum capacity 
and provides lower bounds on the individual-user capac- 
ities. Then we quantify the capacity region that is lost 
when heterodyne or homodyne detection is employed — 
in lieu of the optimum joint measurement — with single- 
use coherent-state encoding. Finally, we derive upper 
bounds on the individual-user capacities, and show that 
they can be achieved — in the limit of high input pho- 
ton numbers — by means of squeezed-state encoding and 
homodyne detection. 



We will begin with the single- mode optical MAC, 
shown in Fig. ^ hi which two senders, Alice and Bob, 
transmit classical information to a common receiver, 
Charlie, by accessing different input ports of a lossless 
beam splitter with transmissivity 77, where < r; < 1. 
The input-output relation for the electromagnetic modes 
associated with this channel is c = — r\ b, where 

a and b are the annihilation operators of Alice's and Bob's 
input modes, and c is the annihilation operator of the 
mode that Charlie measures. In this section, we derive 
the capacity of the optical MAC when Alice and Bob en- 
code complex- valued input messages a and (3 as coherent 
states I a) a <£> \0) b with independent probability densities 
Pa(ck) and pb((3). This encoding puts the c mode in the 
coherent state \^/rja + \/l — r) /3)c, so we will refer to 
this system as the two-user, single-mode, coherent-state 
MAC. 
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A. Quantum MAC Capacity Theorem 



B. Coherent-State MAC Capacity 



The capacity region of a two-user multiple-access chan- 
nel is defined to be the closure of all rate pairs (Ri,R 2 ) 
for which arbitrarily small error probabilities are achiev- 
able in the limit of long codewords 0. Winter's quantum 
MAC capacity theorem |7j gives the capacity region of a 
quantum MAC optimized over arbitrary receiver mea- 
surements, and over codewords that are not entangled 
over multiple channel uses. Winter's result presumes a 
finite-dimensional state space, whereas the Bosonic MAC 
has an infinite-dimensional state space. Nevertheless, we 
shall rely on his result, which can be extended to Bosonic 
channels by means of a limiting argument, see Appendix. 
Thus, the capacity region of the two-user, single-mode, 
optical MAC from Fig. 2] will be taken to be the convex 
closure of all rate pairs (i2i,-Ra) that satisfy the following 
inequalities: 



Ri < j Pb{P)S{p%) dp- J J PA {a) PB {f3)S{p{a, 0)) da d/3, 

(la) 

i?2 < JpA(a)S(p£)da-JJp A ( a )p B {p)S(p(a,P))dadp, 

(lb) 



R1 + R2 < S(p)- II p A (a)p B (0)S(p(a,P))dadp, (lc) 

for some product distribution, p A (a)p B (/3), on Alice's 
and Bob's complex-valued inputs. In these expressions, 
S(-) is the von Neumann entropy, and the average density 
operators are 

Pp = JpA{a)p(a, (3) da, (2) 

Pi = J PB (P)p(a 1 P)dp, (3) 

p= [ fp A (a)p B (P)p(a,P) dad/3, (4) 



where p(a,f3) is the received state given that messages 
a and f3 have been transmitted. The capacity region 
will diverge unless the inputs are constrained, so here we 
assume that Alice and Bob are subject to the average 
photon- number constraints (a^d) < n A and (b^b) < n B , 
respectively. 

Equations constitute the multiple-access version 
of the Holevo-Schumacher- Westmoreland theorem, which 
gives the classical capacity of a single- user quantum chan- 
nel H, U H3- For the coherent-state MAC, p(a,/3) = 
\^Jr\a + a/1 — r](3}(^/r]a + y/1 — 77 (3\ is a pure state, so 
that the second terms on the right-hand sides of these 
equations vanish, and the average density operators in 
the first terms are found by performing the indicated in- 
tegrations in Eqs. 



Suppose that Charlie uses homodyne or heterodyne de- 
tection. These are single-use measurements that may not 
achieve the capacity region of the coherent-state MAC, 
but they are easily realized with existing technology and 
their capacity regions are simple to derive. In particu- 
lar, coherent-state MACs that use homodyne or hetero- 
dyne detection reduce to classical additive Gaussian noise 
MACs: a scalar Gaussian MAC, with noise variance 1/4, 
for homodyne detection, and a 2D white-Gaussian noise 
MAC, with noise variance 1/2 per dimension, for hetero- 
dyne detection. It follows that the capacity region for 
the coherent-state MAC with homodyne detection is the 
set of rate pairs (Ri, R2) that satisfy 6] 



Ri < ~log(l + 4»pu) 

R2 < -log(l + 4(l-77)n B ) 



(5a) 
(5b) 



Ri + R2 < g log(l + 477n A + 4(1 - n)n B ), (5c) 

and the capacity region for the coherent-state MAC with 
heterodyne detection is 



Ri < log(l 
R 2 < log(l 
Ri+R 2 < log(l 



i]n A ) (6a) 
(1 - v)n B ) (6b) 
r\n A + (1 - rf)n B ), (6c) 



when Alice and Bob are subject to the average photon- 
number constraints n A and n B , respectively. 

The preceding two-user results extend easily to the m- 
user coherent-state MAC that employs homodyne or het- 
erodyne detection. Here, the ith transmitter sends coher- 
ent state |ci!j), for 1 < i < m, resulting in the channel's 
output mode being in the coherent state | 2»=i \f^]i (X i)i 
where the transmissivities \r\i\ sum to one. The m-user 
capacity region with homodyne detection is then the set 
of rates (Ri, ■ ■ ■ , R m ) that satisfy the inequalities 



^ J R. 1 <ilog(l + 4^r7. i 



(7) 



for all subsets S C {1, . . . , m}, where (ajflj) < rij is the 
average photon-number constraint on the ith user. Sim- 
ilarly, the capacity region with heterodyne detection is 
given by the inequalities 



(8) 



for all subsets S C {1, . . . , m}. 

The homodyne and heterodyne detection capacity re- 
gions for the coherent-state MAC provide inner bounds 
on that channel's capacity region if no constraints are 
placed on the receiver measurements. We will now 
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FIG. 2: (Color online) Coherent-state capacity region for the 
optical MAC. The capacity region with optimum reception 
(solid line) is given by inequalities 11 U . The capacity regions 
with homodyne detection and heterodyne detection are also 
shown. This figure assumes r\ — 1/2, Ha = 10, and fis = 8. 
Rates are measured in nats, i.e., logarithms are taken base e. 

find the capacity region for the coherent-state MAC — 
without restricting the choice of receiver structure — from 
the previously stated capacity theorem. The channel 
outputs of the coherent-state MAC are the pure states 

p(a,(3) = l-^/rya + Vl — r nP){y/ 7 n a + \A ~ V P\- I* is then 
easy to show that the circularly-symmetric Gaussian dis- 
tributions 

p A (a) = ^ezp(-^), (9) 

ITUA \ UA J 

PBC3) = -^exp(-^V (10) 

are the optimal input distributions, i.e., they maximize 
the right-hand sides of Q for the coherent-state MAC. 
Using these input distributions we find that the capacity 
region of the coherent-state MAC, with optimal (joint 
measurements over entire codewords) reception is the set 
of all rate pairs satisfying 

Ri < 9(wa) (Ha) 
i? 2 < 5((1 - v)n B ) (Hb) 
Ri+R2<g(vnA + {l-v)n B ), (He) 

where g(x) = (x + 1) log(x + 1) — xlog(cc) is the Shan- 
non entropy of the Bose-Einstein probability distribution. 
Figure [21 compares the capacity regions for the coherent- 
state MAC when homodyne detection, heterodyne detec- 
tion, and optimal reception are used. 

Users transmitting information over the optical MAC 
must each contend with interference created by the other 
users who are attempting to access the channel. This 
type of noise, called multiple-access interference, is re- 
sponsible for the pentagonal shape of the capacity re- 
gion seen in Fig. [21 In general, users communicating 
over a multiple-access channel will encounter channel 
noise in addition to multiple-access interference. A two- 



user, single-mode, optical MAC that introduces addi- 
tional white-Gaussian noise can be modeled by the evo- 
lution equation c = ^/rj a + yT — r\ b + £, where £ 
is additive classical zero-mean, complex-valued, white- 
Gaussian noise with variance (|£| 2 ) = N. Our deriva- 
tion of the two-user capacity region for the coherent- 
state MAC generalizes to include the presence of addi- 
tive white-Gaussian noise, with the following result for 
the capacity region: 

Ri < g(vn A + N) - g (N) (12a) 
R 2 < g((l ~ v)n B + N) ~ g(N) (12b) 
Ri+R2< g(vn N + (1 - r))n B + N) - g(N). (12c) 

C. Wideband Capacity 

Now let us turn our attention to the wideband 
coherent-state MAC, in which Alice and Bob may em- 
ploy photons of any frequency, subject to constraints, 
P A and Pb, on their average transmitted powers. For 
a frequency- multiplexed scheme, in which the radian- 
frequency domain is divided into bins of width A = 
2n/T, the channel output for the ith mode is 

Ci = T/fja,i+ y/l-r) h, (13) 

where fij and hi are the input modes at frequency Wi — 
i2ir/T, for i = 1,2,3,..., and the transmissivity 77 is 
frequency independent. The average power constraints 
on Alice and Bob are given by 

Y,^ l (\a l \ 2 )A/2ir < P A (14) 

i 

^2nu l (\p t \ 2 )A/2n<P B , (15) 

i 

where Alice and Bob allocate average photon numbers 
n A (uji) = (\ai\ 2 ) and n B (u>i) = (|ft| 2 ), respectively, to 
frequency-bin u>i. 

We first derive the capacity region of the wideband 
coherent-state MAC with homodyne detection. When 
homodyne detection is employed, the wideband coherent- 
state MAC is equivalent to a set of parallel classical 
MACs with independent zero-mean Gaussian noise, for 
which we have derived upper bounds on the individual 
rates, Ri and i?2, and on the sum rate, Ri + R2, from 
separate Lagrange multiplier calculations. In the limit 
A — > 0, these upper bounds become 

«, + fl2 <y^±£p>^, (16c) 
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with 



r]n A {oj) 



to 



1 n V P A 



for 10 < A^nijPA/h, 



n \ - / \ 1 / tt(1 - y)P B 1 



4' 



for uj < - r/)P B /h, and 



OJ V 



1 /7T[7 ? P A + (1 - 7?)^] 1 

h 4' 



(17) 



(18) 



(19) 



for cj < 4 v /7t[t7P j 4 + (1 — rfyPsl/h, where n^ B 
ryrUM + (1 - i])n B (uj). 

We see that the optimal mean photon number allo- 
cations, n A (u>) and rig (a;), are given by water-filling 
formulas, as is found in classical information theory. 
The rates (|16|) define a pentagonal region which serves 
as an outer bound for the capacity of the wideband 
coherent-state MAC with homodyne detection. To prove 
that this outer bound is, in fact, the capacity region, 
we must show that Eqs. l(T7 )l -lfT9 ^l can be satisfied si- 
multaneously. The average photon number allocations 
(fi A (w), (fi' AB (uj) — T)n A ((j))/(l — rj)) for Alice and Bob, 
achieve the lower-right corner point 



tjPa 



T)P A + {l-ri)P B Kl-v)Pi 



(20) 



of the outer bound. Similarly, ((n^ B (w) — (1 — 
J7)ns(o;))/?7,nB(w)) achieves the upper-left corner. 
Thus, the entire region is achievable and hence is equal 
to the capacity region. A similar derivation shows that 
the wideband coherent-state MAC with heterodyne de- 
tection has a capacity region that is identical to that of 
homodyne detection. 

The preceding two-user wideband results readily ex- 
tend to the m-user wideband coherent-state MAC. Sup- 
pose that the fcth user sends coherent states {lajt^)} 
across the frequency bins {uji}. The channel output for 
the ith-frequency mode will then be the coherent state 
I Sfc=i \/Vk a k,i), where the frequency-independent trans- 
missivities rjk sum to one. The input power constraint on 
the fcth user is 



^/n^(|a M | 2 )A/2^ < P k 



(21) 



for 1 < fc < m. If the receiver uses homodyne or het- 
erodyne detection, then the wideband capacity region is 
defined by the inequalities 



fees V kes 



(22) 



for all5C{l,..., m }. 



In deriving the wideband capacity region for the two- 
user, coherent-state MAC with homodyne or heterodyne 
detection, we first obtained upper bounds on the indi- 
vidual rates R\, i?2, and the sum rate R\ + i?2, and 
then showed that these bounds could be achieved simul- 
taneously. Applying this same procedure to the two- 
user, coherent-state MAC without constraining its re- 
ceiver structure, we have obtained the following capacity 
region, 



Ri < 



R 2 < 



3ft 



7T(1 - 7?)P B 



3h 



Ri + R 2 < 



tt^Pa + (1 - rj)P B ] 



and optimal average photon number allocations, 

1 



r]n A (uj) = 
(1 - T])n B (u) = 



expf ^nhu; 2 /l2r)P A ) - 1 
1 

exp( x /7rftw 2 /12(l - r])P B ) - 1 



(23a) 
(23b) 
(23c) 

(24) 
(25) 



l AB 



M = 



exp y/ntiw*/12[T)P A + (l-T))P B ]) - 1 



(26) 

Equations (|16fl and l|23(l show that optimal reception in- 
creases both the individual rates and the sum rate by a 
factor of 7r/ Vo as compared to what is achievable with 
homodyne or heterodyne detection. The m-user capac- 
ity region for the coherent-state MAC is specified by the 
inequalities 



fees 



< 



fees 



3h 



(27) 



for all S C {1, . . . , m}; once again there is an improve- 
ment factor of 7r/\/3 as compared to homodyne or het- 
erodyne detection. 



III. GAUSSIAN MAC 

Now let us return to the single-mode case and relax our 
assumption that the transmitters use coherent-state en- 
codings, i.e., we will allow them to use non-classical states 
in their quest for the largest possible capacity region. As 
a step toward finding the ultimate capacity region of the 
optical MAC, let us allow Alice and Bob to employ arbi- 
trary Gaussian states, instead of just coherent states. 



A. Holevo-Sohma-Hirota MAC 

It is useful to begin by deriving the capacity re gion 
for the two-user version of the Holevo-Sohma-Hirota fijl 
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(HSH) channel model. Consider a mode with annihila- 
tion operator a = a\ + ia 2 that is in a zero-mean, Gaus- 
sian state, /5(0), with quadrature-component covariance 
matrix 



V 



(aia 2 + a 2 ai)/2 



[a\a 2 + a 2 ai)/2 



We define a multiple access channel model in which Alice 
and Bob send classical messages a and /?, subject to input 
constraints 



\a\ 2 PA {a)da = N A , (29) 
\p\ 2 p B (/3)d/3 = N B , (30) 



and Charlie receives the shifted version of the initial 
state, viz., p{a,/3) = D{a + p)p(0)Di(a + 0), where 
l)(7) = exp(7a^ — 7*a) is the displacement operator. 
From the quantum MAC capacity theorem Q, the ca- 
pacity region of the two-user HSH MAC is given by the 
convex hull of all rate pairs (Ri, R 2 ) satisfying 

Ri < S(p A ) - S(p(0)) (31a) 
R 2 < S(p B ) - S(p(0)) (31b) 
Ri+R2<S(p AB )-S(p(0)), (31c) 

for some product distribution p A (a)p B ((3), where the av- 
erage density operators are 

PA= [p A {a)D(a)p(0)&(a)da, (32) 
Pb = JpB((3)D(P)p(0)DHP)d(3, (33) 
Pab = [p A (a)p B {P)p{a,P) dad/3. (34) 



To evaluate this capacity region, we will first maximize 
the rate upper bounds for Ri, R 2 , and R\+R 2 separately. 
Then we will then show that the region described by these 
maximum rates is achievable. 

To maximize the R\ upper bound in l|31a|) . we follow 
the proof of the HSH capacity theorem [llj; the same 
derivation will also apply to the R 2 upper bound. For any 
input distribution p A {a) that satisfies constraint l|2*9"|l . let 
p A (a) be the zero-mean Gaussian distribution with the 
same second moments as p A (a). Then, p A {a) satisfies 
constraint (|29|l and 



PA 



PA 



(a)D{a)p(0)D\a) da 



(35) 



is a Gaussian state. If F(a, a') is any second-order poly- 



nomial in (a, a'), then 

tr[p A F(a, a f )] = J p A {a) tr[D(a)p(0)D^ (a)F(a, a r )} da 

(36) 

p A (a) tr[p(0)F(a + a, + a*)) da 

(37) 

p A (a) tr[p(0)F(a + a, + a*)] da 

(38) 

= tr[p A F{a > a^)}. (39) 

Thus, p A and p A have the same second moments, and it 
follows that S(p A ) > S(p A ), i.e., we can restrict our at- 
tention to Gaussian input distributions in trying to max- 
imize the i?i upper bound. 

When the input distribution p A (a) is Gaussian, the 
rate upper bound for R\ can be expressed as 

S(p A )-S(p(0)) = g(2\V+V a \ 1 ^-l/2)-g{2\V\ 1 / 3 -l/2), 

(40) 

where the quadrature-component covariance matrix of 
p A {a) is 



V a = 



v& v 2 a 



(41) 



Thus, the optimization problem we need to solve is 
maxy Q |V + V a \, subject to the positive semidefinite and 
input power constraints 



V a >0, 
ix(V a ) = V? + V 2 a = N A . 



(42) 
(43) 



This constraint region is the interior of a circle in the 
Vf — plane, which has the following polar-coordinate 
parameterization, 

„, N A 



2 ' 



V,%=rsm0, V 2 a 



where < r < N A /2 and < 9 < 2ir. Now, write 

\v + v a \ = (Vi + vf)(y 2 + v 2 a ) - (v 12 + v? 2 f 

2 



V 1 + V 2 + N A 



{Vr. 



Na 

2 ' 
(44) 

(45) 

2 



+ r cos ( 



(46) 



In terms of r and 8, our maximization problem then be- 
comes 



maxlV + V a \ 



V x + V 2 + N A 



mm 

rfi 



2 

V 2 



Vi 



-V r . 



(47) 
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This maximization has two different solutions, de- 
pending on whether or not the point ((Va — Vi)/2, 
— V12) lies in the radius- N A /2 circle whose center is at 
the origin. If ((V 2 — Vi)/2, — Via) lies in this circle, then 
the minimum on the right-hand side of (|47() is zero. If 



((V 2 — Vi)/2, — V12) lies outside this circle, then a simple 
geometric calculation gives the minimum on the right- 
hand side of 147(1 . We thus obtain the maximum individ- 
ual rates 



-Rmaxi = maxS(pA) - S(p(0)) 



{ g{V l + V 2 +N A -\)-g{2\V\ 1 / 2 ~\) 1 

for N A >((V 1 -V 2 ) 2 + 4V 2 2 ) 1 / 2 , 



9\ 2 



[(Vi +V 2 + N A )/2] 2 - - V 2 )/2] 2 + W2 - Na/2)' 

for < ((Vi - V 2 ) 2 + 4V 1 2 2 ) 1 /2 

I 



1/2 



(48) 



(49) 



-0(2|V| 



1/2 _ i\ 
21 ' 



A similar expression holds for the maximum rate i? ma x2 • 
The fact that capacity is given by two different expres- 
sions depending on whether input constraints satisfy a 
certain inequality is referred to as a noncommutative gen- 
eralization of waterfilling in [TT| . 

To maximize the sum-rate upper bound, we follow the 
same approach. It is again sufficient to consider Gaussian 
input distributions p A (ct) and Pb{P), so our maximiza- 
tion problem is maxy ,Vfl \V + V a + Vp\, subject to the 
positive semidefinite and input power constraints 



v a > 0, Vp> 0, 



(50) 



tr(V tt ) = V X Q + V 2 a = N A , tr(V» = Vf + Vg = N B . 

(51) 

This constraint region is the interior of two circles whose 
polar-coordinate parameterizations are 



Va 



V a = 



r A cose A r A sm9 A \ Na 
r A smv A —r A cos6 A J 2 



r b cos 6b tb sm " B 
rssmOB —tbcosOb 



Nn 

-fl, (53) 



for < r A < N A /2, < r B < N B /2, < 6 A < 2vr, and 
< Ob < 2ir, with / being the 2x2 identity matrix. 
This parameterization allows us to write 

\V + V a + V p \ = (V 1 + V? + Vf)(y 2 + V? + Vf) 

-(V^ + VZ + vfz) 2 (54) 
V1 + V2 + N A + N B 



V x -V 2 



r A cos v A + tb cos u b 



-(V 12 +r A sme A +r B sm9 B ) 2 - (55) 



Thus, we have 



max IV + V a + Vb\ 

VcV/3 



Vi + V 2 + N A + N B 



— mm 

r A ,rB,0A- 



V2-V1 



r A cos o A — tb cos Vb 



+(-V*i2 - r A smQ A -r B sin 0b) 2 



(56) 



The second term on the right in (|56J) is the minimum 
squared distance between the points ((V 2 — Vi)/2, —Via) 
and (r A cosd A + vbCOsUb, r A sin6 A + t\b siring). If 
((Va - Vi)/2,-Via) lies in the radius-^ + N B )/2 cir- 
cle that is centered at the origin, then the second term 
vanishes. Otherwise, a simple calculation gives this min- 
imum distance. As a result we find that the maximum 
sum rate is, 



7 



-R max i2 = max S(pab) - S(p(0)) 

V a ,Vp 



(57) 



f g{Vi + V 2 + N A + N B - |) - g(2\V\ 1 / 2 - \) , 

for N A + N B > ((Vi - V 2 ) 2 + W 2 ^ 1 ' 2 



9\ 2 



[(Vi + V 2 + N A + N B )/2] 2 - [y[(Vi - V 2 )/2} 2 + V 2 2 - (N A + N B )/2 

-g(2\V\^ 2 -\), 

for N A + N B < ((Vi - V 2 ) 2 + W^) 1 ' 2 . 



1/2 



(58) 



We claim that the two-user capacity region for the HSH 
MAC with initial-state quadrature-component variance 
matrix V and input constraints N A and N B , is the region 
defined by the inequalities 

R\ < -Rmaxl, R 2 < -Rmax2, and R\ 2 < i?maxl2- 

(59) 

To verify this claim, we will show that the corners of this 
region are achievable. The capacity region result then 
follows by time-sharing. 

Let the point ((V 2 — Vi)/2,—Vi 2 ) have coordinates 
(rv,6V) and suppose that N B > N A . To show that 
the lower corner (i? max i, i? ma xi2 - R ma , x2 ) is achievable, 
we need to find points (r A ,6 A ) and (r B ,9 B ) that si- 
multaneously minimize the distance between (ry,9v) 
and (r A ,9 A ) and the distance between (rv,&v) and 
(t A ,0a) + (r B ,9 B ). Similarly, to show that the upper 
corner (i? max i2 - Rmaxi, Rmax2) is achievable, we need to 
minimize the distance between (ry, Ov) and (r B , 9 B ) and 
the distance between (ry,0v) and (r A ,9 A ) + (r B ,9 B ). 
There are four cases to consider: see Fig. [3] For each 
case, we list the coordinates (r A ,9 A ) and (r B ,9 B ) corre- 
sponding to the capacity-achieving input distributions. 

• Case I. 

— lower corner: (r A ,9 A ) — (ry,9v) and 
(r B ,0 fl ) = (O,O) 

— upper corner: (r A , 9 A ) = (0, 0) and (r B , 9 B ) = 
(rv,0 v ) 

• Case II. 

— lower corner: (r A ,9 A ) = (N A /2,9y) and 
(r Bl 9 B ) = (r v -N A /2,0 v ) 

— upper corner: (r A , 9 A ) — (0, 0) and (r B , 9 B ) = 
(rv,0v) 

• Case III. 

— lower corner: (r Al 9 A ) — (N A /2,9y) and 
(r Bl 9 B ) = (rv-N A /2,6 v ) 

— upper corner: (r A ,9 A ) — (ry — N B /2, 6y) and 
(r B ,9 B ) = (N B /2,9 v ) 




FIG. 3: (Color online) Regions used to complete the HSH 
MAC capacity-region proof. Case I: rv < N A /2. Case II: 
N A /2 <r v < N B /2. Case III: N B /2 < r v < (N A + N B )/2. 
Case IV: r v > (N A + N B )/2. 

• Case IV. 

— lower corner: (r A ,9 A ) = (N A /2,9y) and 
(r B ,0 B ) = (N B /2,9 V ) 

— upper corner: (r A ,9 A ) — (N A /2,9y) and 
(r B ,9 B ) = (N B /2,9 v ). 



B. Gaussian MAC Capacity 

We now apply the capacity result derived in the previ- 
ous section to the two-user optical MAC in Fig. Alice 
and Bob encode their classical messages a and (3 using 
input states of the form 

PA (a) = b{a)p A {Q)b\a) (60) 
PB ((3) = D(p)p B (0)DHp), (61) 

where p A (0) and p B (0) are zero-mean Gaussian states 
with quadrature-component covariance matrices V A and 
V B , respectively. This is a modulation code for which the 
coherent-state encoding is the special case in which p A (0) 
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and pb(0) are vacuum states. Charlie receives the output 
ensemble {pa(o)Pb(P), £{Pa(&) <8> Pb (/?))}> where the 
channel output £(pa(cx) ® pb{P)) is the Gaussian state 
with mean ^/rj a + y/1 — r\ f3 and covariance matrix rj Va + 
(1 — T))Vb- The capacity of this Gaussian MAC, with 
input mean photon number constraints ua and n^, is 
the HSH capacity region found in the previous section 
with 

V = V V a + (1-v)Vb, (62) 
N A = v(n A - Vx - V 2 A + 1/2) , (63) 
N B = {1- rf) (n B - V? - V 2 B + 1/2) . (64) 

For ua and Ub sufficiently large, this capacity region is 
the set of rate pairs that satisfy 

Ri < 9{vn A + (1 -rj)(Vf + V 2 B - 1/2)) 

-g(2\V\^ 2 -l/2) 7 (65a) 

i? 2 < g(v (Vf + V 2 A - 1/2) + (1 - n)n B ) 

-g{2\V\ l ' 2 -l/2), (65b) 

R1 + R2 < 5 (?7n J 4 + (l-^)n B )-.g(2|l/| 1 / 2 - 1/2) . (65c) 

When Va = Vb = 1/4, Alice's and Bob's initial states 
are vacuum states, hence they are employing coherent- 
state encoding and Eqs. I|65ll reduce to the coherent-state 
formulas in Qllfl. As shown in Fig. EI it is possible to find 
Va and Vb, for example, 



V A = V B 



1/32 
2 



(66) 



when 77 = 1/2, ua — 10, and fis = 8, such that the Gaus- 
sian MAC capacity region is larger than the coherent- 
state MAC region. Numerical search over the space of 
possible covariance matrices is one way to further enlarge 
the capacity region beyond that achieved by this exam- 
ple. In the next section, we derive a result that implies 
the maximum individual rates achievable over the Gaus- 
sian MAC when Alice, say, is allowed to choose the opti- 
mal input covariance matrix Va corresponding to Bob's 
covariance matrix Vb- In Section^] we show that trans- 
mitting Gaussian states is asymptotically optimal in the 
limit of large ua and Tib ■ 



IV. ANISOTROPIC GAUSSIAN-NOISE 
CAPACITY 

In this section, we generalize previous work [B S El- 13 
on single-user lossy Bosonic channels with Gaussian ex- 
cess noise to include anisotropic (colored) noise. In this 
section, the channel model we shall consider is the trace- 
preserving completely-positive (TPCP) map, £Y b {-), as- 
sociated with the evolution from input mode a to output 



mode c = ^/rj a + \J\ — 77 b, when the noise mode, 6, is in 
a zero-mean Gaussian state, pb, with quadrature covari- 
ance matrix Vb- Let n& denote the mean photon number 
of the Gaussian noise state pb- 

In seeking the capacity of this channel, we shall assume 
that the conjecture about the minimum output entropy 
of the thermal-noise (isotropic-Gaussian noise) channel 
Q is correct. This conjecture states that the minimum 
output entropy of the thermal-noise channel £ (•), which 
has Vb = (2ht + 1H/4, where / is the 2x2 identity 
matrix, is given by 



.S(£(p))=g((l-r 1 )n T ) 



(67) 



Presuming the correctness of this conjecture, we now 
have the following theorem. 

Theorem 1 The classical capacity of the Gaussian-noise 
channel £^ b is given by 

C = g( V n + (1 - rj)n b ) - g ((1 - rj) (2\V b \ 1/2 - 1/2)) , 

(68) 

for input mean photon numbers n > fithresh; where 



™thr 



sh = - ((v; - vif + 4Vi'|) 1/2 + v + V 2 - l -, 



V = rjV + (1 - v)v b , 

1 / |/J + H 2 2Im(/w) 



V 



4 l2Im(/^i/) \fi 



(69) 
(70) 

(71) 



and the parameters p and v are chosen such that the 
squeeze operator S(z) whitens the Gaussian state pb- 

For sufficiently large input mean photon number n, l|68|) 
gives the classical capacity of the single-user Gaussian- 
noise channel. 

Proof We begin by establishing an upper bound on 
the capacity. By the HSW theorem, 



C < max S $X 6 (Pi) " min V Pj S(^(ft)) 
!/■ ,P • \ / {Pj >Pj } j 



< g(j]n + (1 - rj)n b ) - min S(£^ b (£ .)). 



(72) 
(73) 



As sketched in Fig. 01 we can use the unitary squeeze op- 
erator S(z) to find a thermal-noise channel, with TPCP 
map £(■), whose output minimum output entropy is equal 
to that of our anisotropic noise channel. [In essence, 
this is the quantum equivalent of the noise-whitening ap- 
proach to communication through colored noise that is 
employed in classical communication theory.] The aver- 
age noise-photon number, fix, of this equivalent channel 
is 



fi T = 2\V b \ 



1/2 



1/2, 



(74) 
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FIG. 4: Equivalent thermal-noise channel £% T from a to 
c'. The input mode a' is in state p' , and the noise op- 
erator b is in a thermal state with mean photon number 
fiT = 2jVi,| 1 ^ 2 — 1/2. The original Gaussian-noise chan- 
nel takes input a to output c. For the squeeze operator 



S(z) = exp[(z*a 



)/2], we use the parameterization 



/i = coshr and v = e lS sinhr, where z = re 19 . Squeezed 
vacuum states are defined as |0; z) = S(z) 0). 



g((l — 7])n,T), in accord with the capacity conjecture for 
the thermal-noise channel £| . A more interesting special 
case occurs when pi, = \0;z)(0;z\ is a squeezed state, 
with > 0, i.e., a pure-state anisotropic Gaussian noise. 
Here we find 



V' = V = v b 



\fi-v\ 



—2 Im(/ij/) 



-2Im(/i^) \fi 



which yields 



(78) 



(79) 



for n > Uthresh = ImH/^ + W\ ■ Note that this ca- 
pacity is higher than that of the thermal-noise channel 
with the same value. In other words, phase-sensitive, 
pure-state Gaussian noise enhances, rather than degrades 
channel capacity for n > nthresh- 



CAPACITY OUTER BOUND 



which, when used in conjunction with i|73[l and our min- 
imum output entropy conjecture, shows that the right- 
hand side of (|68|l is an upper bound on the channel ca- 
pacity. 

To show that the right-hand side of is also a 

lower bound on the channel capacity when fi > nthresh, 
we evaluate the information rate achieved by a single- 
use squeezed-state code. Let fP a = |0;— z)(0;— z\ be the 
zero-mean squeezed state whose quadrature-component 
covariance matrix is given by (|71(l . Consider that ran- 
dom code in which we transmit the displaced squeezed 
states, 

p a (a) = D{a)p a b\a), (75) 

that are selected with a zero-mean Gaussian probabil- 
ity density function whose quadrature-component covari- 
ance matrix is denoted V a . Imposing the average photon 
number constraint, (a^a) < n, assumin g th at n > nthresh; 
and applying the HSH capacity result |TTI| - we find that 
there is a squeezed-state code whose information rate 
equals the right-hand side of (36). This implies that 

C > g(r,n + (1 - rj)n b ) - g ((1 - r,) (2|H| 1/2 - l/a)) . 

(76) 

Equations (41) and (44) provide coincident upper and 
lower bounds on the capacity, when n > n t hresh, hence 
the proof is complete. 

There are two special cases of this theorem that are 
worth discussing. First, it is easy to see that when 



£^ b reduces to the thermal-noise channel £. Theorem 1 
then predicts nthresh = and C = g{rjfi + (1 — rfjur) — 



Achieving the ultimate capacity region of the optical 
MAC may require the use of non-Gaussian states, so 
the capacity of the Gaussian MAC is still only an in- 
ner bound on this region. In this section, we develop an 
outer bound on the ultimate capacity region of the opti- 
cal MAC. Let Alice and Bob use input states — averaged 
over their respective random-coding ensembles — pA and 
pB that are subject to the average photon number con- 
straints fiA and fiB- Because von Neumann entropy is 
invariant to mean fields, we know that the optimum pA 
and pB will be zero-mean-field states. This, in turn, im- 
plies that (c^c) = rjhA + (1 — i])^B^ from which it is easily 
shown that 

Ri + R 2 < S{£{p A ® Pb)) < g (wa + (1 - r))n B ) ■ (80) 

The sum-rate upper bound in l|8UI) coincides with the 
coherent-state MAC result appearing in jllcf . Hence, we 
have shown that the sum rate for the capacity region is 
achieved by coherent-state encoding in conjunction with 
optimum (joint-measurement) reception. More gener- 
ally, the Gaussian-state encoding is a sum-rate-achieving 
code in the above-threshold regime, i.e., H65cfl coincides 
with (|%U|l . whenever £(pa(0) ® Pb(0)) is pure. Moreover, 
from (|6cll it can be shown that heterodyne reception is 
asymptotically optimum for the sum rate in the limit 
r\n A + (1 - r))n B — > 00. 

To upper bound the individual rates R\ and R2, con- 
sider a super receiver that has access to both output 
ports of the beam splitter representing the optical MAC. 
This super receiver can invert the unitary beam splitter 
transformation to undo the effects of the optical MAC. 
Thus, the individual rate upper bounds reduce to single- 
user Holevo informations, and we have the upper bounds 
Ri < s(nyi) and R 2 < g{fiB)- Our optical MAC results 
are illustrated in Fig. \E\ Here we have plotted the sum 
rate for a two-user, single-mode, quantum optical MAC 
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ult. cap. outer bound 

Gaussian encoding 

coh. -state cap. 

— ■ ■ coh. -state +het. 




FIG. 5: (Color online) Ultimate capacity region of the two- 
user, single-mode optical MAC. Inner bounds from coherent- 
state and Gaussian-state encodings, and outer bounds given 
by Ri < g{n A ), R2 < g{n B ), and i?i + R 2 < g(j]n A + (1 - 
rj)fiB) are shown. The Gaussian-state capacity region is eval- 
uated with input variance matrices Va and Vb given by 1661 . 
This figure assumes 77 = 1/2, ua = 10, and ub = 8. Rates 
are measured in nats, i.e., logarithms are taken base e. 



with rj = 1/2, n A = 10, and jig = 8, along with the 
capacity region for heterodyne detection, the individual 
rate limits for coherent-state encoding, and the individual 
rate limits for the Gaussian-state encoding from Eq. (|66H • 
We have presented codes which achieve the sum-rate 
upper bound, but it is unknown exactly how far we can 
reach into the corners of the outer bound region. One 
thing we can demonstrate is that the individual rate up- 
per bounds are asymptotically achievable in the limit of 
large n A and ub- Let Alice transmit real- valued classi- 
cal messages a\ using squeezed states |ai;z) excited in 
the first quadrature with squeeze parameter z > 0. Let 
Bob transmit the zero- mean squeezed state |0; Z) with 
squeeze parameter Z = sinh _1 (y 'Ub), i-e., Bob squeezes 
as hard as possible, given his average photon number 
constraint. A rate of 



Ri 



4(n A 



sinh 2 z) 



+ (1 - r))e- 2Z /r) 



(81) 



is achieved if Charlie uses homodyne detection to decode 
Alice's message. After substituting in Alice's optimal 
value for her squeeze parameter, z = log(2n>i + l)/2, 
and performing several applications of L'Hopital's rule, 
we obtain the ratio 

lim lim 1 

ua^oo n B ^<x> g[n A j 

\ log (1 + Ae 2z {n A - sinh 2 z)) 



lim 
lim 

UA — >OC 
1. 



g(n A ) 



log(l + 2n A ) 



(82) 

(83) 
(84) 



Thus, this squeezed-state code with homodyne detection 
is asymptotically optimal for large input photon numbers 
n A and t%b- For the special case 77 = 1 , Bob is irrele- 
vant, and the above argument shows that the squeezed- 
state/homodyne code is asymptotically optimal for the 
single-user lossless Bosonic channel. 



VI. CONCLUSIONS 

We have derived the capacity region of the Bosonic 
multiple-access channel that uses coherent-state encod- 
ing. Single-mode and wideband transmitters were con- 
sidered, and in both cases optimum (joint measurements 
over entire codewords) reception was shown to outper- 
form receivers that employed homodyne or heterodyne 
detection. Coherent-state encoding with optimum recep- 
tion was shown to achieve the sum-rate bound on the 
ultimate capacity region of the optical MAC. In the limit 
of high average photon numbers, the ultimate single-user 
rates can be achieved with squeezed-state encoding and 
homodyne detection. 
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APPENDIX 

We apply the capacity theorem derived in based on 
the following argument. Suppose that the transmitters 
used by Alice and Bob employ states containing no more 
than K photons, where 



K > max{l, n A + ub} ■ 



(A.l) 



As their states may be indexed by complex-valued pa- 
rameters a and f3, over which we can do random cod- 
ing, the result described by QJ specifies the achievable 
rate region within this restricted finite-dimensional state 
space. The right-hand sides of — when maximized 
over product distributions that respect the n A and t%b 
constraints — are monotonically expanding achievable 
rate regions with increasing K . Moreover, the achievable 
rate region for any K is outer bounded by the results 
we derive in Section assuming the full Hilbert space 
is employed. For fixed n A and ub, the impact of the 
truncation to no more than K photons will be negligible, 
under the condition 1A.1|) . and will vanish as K — > 00. 
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